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We present a numerical study of the isotope effect on the angle resolved photoemission spec- 
tra (ARPES) in the undoped cuprates. By the systematic-error-free Diagrammatic Monte Carlo 
method, the Lehman spectral function of a single hole in the tt't" — J model in the regime of inter- 
mediate and strong couplings to optical phonons is calculated for normal and isotope substituted 
systems. We found that the isotope effect is strongly energy-momentum dependent, and is anoma- 
lously enhanced in the intermediate coupling regime while it approaches to that of the localized hole 
model in the strong coupling regime. We predict the strengths of effect as well as the fine details of 
the ARPES lineshape change. Implications to the doped case are also discussed. 
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It is a subject of considerable debates for many years 
whether the electron-phonon interaction (EPI) is essen- 
tial for the physics of undoped and heavily underdoped 
high T c superconductors [1 S H S H H . Early study 
of the small isotope effect (IE) on T c at optimal doping 
together with the no prominent phonon features in the 
temperature dependence of the resistivity p(T) had led 
to the conclusion that the phonons are a by-player in 
high T c cuprates. However, the roles of EPI have been 
recently studied intensively in terms of the neutron scat- 
tering, Raman scatteringjthe IE on T c and the superfluid 
density p s , and ARPES Q. Especially ARPES provides 
the direct information on the single-particle Green func- 
tion, which has revealed the "kink" structure of the elec- 
tron energy dispersion 40-70 meV below the Fermi en- 
ergy. The origin of this kink is naturally interpreted as 
the coupling to some bosonic modes. The magnetic res- 
onance mode and the phonon modes are the two major 
candidates, and the IE on ARPES should be the smoking- 
gun experiment to distinguish between these two. Gweon 
et al. performed the ARPES experiment on O 18 - 
replaced Bi2212 at optimal doping and found an appre- 
ciable IE, which however can not be explained within 
the conventional weak-coupling Migdal-Eliashberg the- 
ory. Namely the change of the spectral function due to 
18 -replacement has been observed at higher energy re- 
gion beyond the phonon energy (~ 60meV). This is in 
sharp contrast to the weak coupling theory prediction, 
i.e., the IE should occur only near the phonon energy. 
Hence the IE in optimal Bi2212 remains still a puzzle. On 
the other hand, the ARPES in undoped materials, i.e., 
the single hole Lehman spectral function (LSF) doped 
into the Mott insulator, y| has recently been understood 
in terms of the small polaron formation [gL [ToL ITT| ■ It has 
been revealed that the energy dispersion of the extended 
t — J model, tt't" — J model [IJ survives as the cen- 
ter of mass position of the broad Franck-Condon peak 
(FCP) of phonon side-bands even at strong EPI, while 
the weight Z of the zero-phonon line and its dispersion 
is very small. Recently, decoupling of chemical potential 



from observed in ARPES resonance, predicted by strong 
EPI scenario, was confirmed experimentally [j| and, thus, 
phonon origin of ARPES broadening seems to be likely. 

In addition to high-T c problem, strong EPI mechanism 
of ARPES spectra broadening was considered as one of 
alternative scenarios for diatomic molecules ^| , colossal 
magnetoresistive manganites jl4j| . quasi-one-dimensional 
Peierls conductors 0], an d Verwey magnetites (l6| . 
Therefore, exact analysis of the IE on ARPES at strong 
EPI is of general interest for conclusive experiments in a 
broad variety of compound classes. 

In this Letter we present a study of the IE on the single 
particle LSF of a hole strongly interacting with phonons 
in the tt't" — J model, which is equivalent to study of 
ARPES in undoped cuprates |4j . As mentioned above, it 
is essential to compare experiment in undoped systems 
with the present DMC results, where theory can offer 
quantitative approximation- free results. We calculate IE 
on LSF by exact DMC method at zero temperature in 
nodal (7r/2,7r/2) and antinodal (ir, 0) points of the Bril- 
louin zone for realistic parameters of the tt't" — J model in 
the intermediate and strong EPI regime. We find that IE 
is anomalously enhanced in the intermediate EPI regime 
but approaches simple analytic estimates in the strong 
EPI regime. We analyze the features of the FCP which 
are most sensitive to IE at different EPI regimes. 

In the standard spin-wave approximation for two- 
dimensional tt't" — J model, which was shown to 
be a good approximation for small exchange integrals 
J/t < 0.4 [13, a hole with dispersion [l^] e(k) = 
4t' cos(k x ) cos(ky) + 2t" (cos(2k x ) + cos(2k y )) propagates 
in the magnon and phonon (annihilation operators are 
Qfk and 6k, respectively) bathes 

Hi j = ]T e(k)hlh k + J2 Wk4«k + O E b l b * (!) 



with magnon dispersion lo^ = 2Ja/1 — 7^, where 7k = 
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(cosfcr + cosfcj,)/2. The hole is scattered by magnons 
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where Mk, q is the standard vertex [17]]. We chose the 
simplest Holstein short-range EPI Hamiltonian 
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where a is the momentum and isotope independent cou- 
pling constant, M is the mass of the vibrating lattice ions, 
and f2 is the frequency of dispersionless phonon (Planck 
constant Ti is set to unity). We introduce dimensionless 
coupling constant A = 7 2 /4tf2 which is, in contrast to the 
standard Holstein constant 7 = aj ^/{2MQ), is invariable 
quantity for the simplest case of IE. Indeed, assuming the 
simplest natural relation Q ~ 1 j\TM between phonon fre- 
quency and mass, we find that A does not depend on the 
isotope factor k- 1so = Q/Qq = y/ Mq/M, which is defined 
as the ratio of phonon frequency in isotope substituted 
n and normal Qq system. 
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FIG. 1: Hole LSFs at k = (n/2,n/2) (solid line) and k = 
(•7T, 0) (dashed line) for different couplings. 
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FIG. 2: Low energy part of LSF for A = 0.62 (solid line), 
A = 0.69 (dashed line), and A = 0.75 (dotted line) at nodal 
(a) and antinodal (b) points. 

We chose adopted parameters of the tt't" — J model 
which reproduce the experimental dispersion of a quasi- 
particle (QP) J/t = 0.4, t'/t = -0.34, and t"/t = 



0.23 . The frequency of the relevant phonon [4j is set to 
fl/t = 0.2 and the isotope factor Kj so = ^/16/18 corre- 
sponds to substitution of O 18 isotope for O 16 . 

We use DMC method [3 which is the only 

available method for approximation-free study of excited 
states of problem llll.'il) for the macroscopic system at 
strong EPI. Non-crossing approximation for phonon vari- 
ables 20] is shown to be invalid approximation for strong 
EPI 9] while exact diagonalization studies of small clus- 
ters, although account for correct treatment of phonons, 
imply a discrete spectrum and, thus, tiny changes of the 
FCP lineshape for k w 1 are not reliable. The only study 
of IE by exact method [2lJ, i.e. path integral quantum 
Monte Carlo algorithm |22|. does not addresses realistic 
tt't"— J model, where both phonon and magnon variables 
have to be summed simultaneously and does not study 
ARPES directly. On the other hand, DMC method al- 
lows to treat ARPES directly . We use recently devel- 
oped Stochastic Optimization method |l8j| which avoids 
regularization and artificial broadening of LSF peaks. To 
sweep aside any doubts of possible instabilities of analytic 
continuation, we calculate the LSF for normal compound 
(n nor = 1), isotope substituted (Ki so = y/16/18) and 

"anti-isotope" substituted (K an t = \/l8/16) compounds. 
Monotonic dependence of LSF on k ensures stability of 
analytic continuation and gives possibility to evaluate the 
error-bars of a quantity A using quantities ^li SO — -4 n0 r, 

«^4nor **4ant ) and (w4i so -4ant)/2. 

Figure presents LSFs in nodal and antinodal points 
for values of EPI which are larger than the critical 
self-trapping (ST) coupling A cr w 0.58, i.e. coupling 
where fast transformation to regime of strong EPI oc- 
curs |2![2l|. Since LSF is sensitive to strengths of EPI 
only for low frequencies we concentrate on the low energy 
part of the spectrum (Fig.|5J. It is seen that near the ST 
crossover the LSF in the nodal point quickly changes with 
coupling while FCP in antinodal point gradually broad- 
ens. The bandwidth of FCP dispersion is estimated as 
Wfcp ~ 0.9. 

Figures El and ^ show IE on the hole LSF for dif- 
ferent couplings in nodal and antinodal points, respec- 
tively. The general trend is a shift of all spectral fea- 
tures to larger energies with increase of the isotope mass 
(k < 1). One can also note that the shift of broad FCP 
is much larger than that of real narrow QP peak. More- 
over, for large couplings A the shift of QP energy ap- 
proaches zero and only decrease of QP spectral weight 
Z is observed for larger isotope mass. On the other 
hand, the shift of FCP is not suppressed for larger cou- 
plings. Except for the LSF in nodal point at A = 0.62 
(Fig. b) , where LSF still has a notable weight of QP 
(5-functional peak, there is one more notable feature of 
the IE. With increase of the isotope mass the height of 
FCP increases. Taking into account the conservation law 
for LSF J^°° Lk(u>) = 1 and insensitivity of high energy 
part of LSF to EPI strength (Fig. 0) , the narrowing of 
the FCP for larger isotope mass can be concluded. To 
understand the trends of the IE in the strong coupling 
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regime we analyze the exactly solvable independent os- 
cillators model (IOM) |25j. More rigorous Lang-Firsov 
transformation is not required since it is not valid for 
adiabatic case Vt/t = 0.2 and in the intermediate cou- 
pling regime pH while in the strong coupling regime, 
as it is shown below, IOM shows quantitative agree- 
ment with DMC data. Indeed, coherent QP bandwidths 
PF QP (A = 0.62) re 3.2 * 10~ 3 , W QP (X = 0.69) «6* 10~ 4 , 
and Wqp(X — 0.75) pa 4 * 10~ 4 are negligibly small in 
comparison with FCP bandwidth Wfcp ~ 0.9. The 
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FIG. 3: Low energy part of hole LSFs in the nodal point 
at different couplings (a, c, e): normal compound (solid line), 
isotope substituted compound (dotted line) and "antiisotope" 
substituted compound (dashed line). Insets (b, d, e) show low 
energy real QP peak. 
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FIG. 4: Low energy part of hole LSFs in the antinodal point. 
See caption of Fig. [3] 



LSF in IOM is the Poisson distribution 

LH = ex P Ko/«] E ^Tp^H 



(4) 



where £o = 7o/^o = 4iA/f2o is dimensionless coupling 
constant for normal system and G k j{lu) = S[lu + 4iA — 
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FIG. 5: Energies of ground state and broad peaks (a) for nor- 
mal (tringles), isotope substituted (circles) and "antiisotope" 
substituted (dimonds) compounds. Comparison of IOM es- 
timates (lines) with DMC data in the nodal (squares) and 
antinodal (diamonds) points: shift of the FCP top (b), FCP 
leading edge at 1/2 of height (c), and FCP leading edge at 
1/3 of height (d). 



QqkI\ is the (5-function. The properties of the Poisson 
distribution quantitatively explain many features of the 
IE on LSF. 

The energy uqp = —At\ of the zero-phonon line I = 
in (0J depends only on isotope independent quantities 
which explains very weak isotope dependence of QP peak 
energy in insets of Figs. 13141 Besides, change of the 
zero-phonon line weig ht obeys relation zj®/Z§& = 
cxp [— £o(l — K )/ K ] m IOM. These IOM estimates agree 
with DMC data within 15% in the nodal point and 
within 25% in the antinodal one. IE on FCP in the 
strong coupling regime follows from the properties of zero 
M = J+™L{u)(L) = 1, first M l = /+~ uL(w)dj = 0, 
and second Mi = f_°° lu 2 L(uj)<Lj — k£o^o moments of 
shifted Poisson distribution J2J ■ Moments M and M 2 es- 
tablish relation V = hf^/h^ = 1/s/k ps 1.03 between 
heights of FCP in normal and substituted compounds. 
DMC data in the antinodal point perfectly agree with the 
above estimate for all couplings. This is consistent with 
the idea that the anti-nodal region remains in the strong 
coupling regime even though the nodal region is in the 
crossover region. In the nodal point DMC data well agree 
with IOM estimate for A = 0.75 {V ps 1.025) whereas at 
A = 0.69 and A = 0.62 influence of the ST point leads 
to anomalous values of V: V ps 1.07 and T> ps 0.98, re- 
spectively. Shift of the low energy edge at half maximum 
Ax/2 must be proportional to change of the root square 
of second moment Ayj^- = y^^oll — V^\- As we found 
in numeric simulations of Eq. (0J with Gaussian func- 



tions [23 Gn.iiy), relation A x / 2 ~ Ay^/2 is accurate 
to 10% for 0.62 < A < 0.75. Also, simulations show 
that the shift of the edge at one third of maximum Ax/3 
obeys relation A 1 / 3 ps A^jj^. DMC data with IOM es- 
timates are in good agreement for strong EPI A = 0.75 
(Fig.lSJ). However, shift of the FCP top A p and A 1 / 2 are 
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considerably enhanced in the self-trapping (ST) transi- 
tion region. The physical reason for enhancement of IE 
in this region is general property regardless of the QP 
dispersion, range of EPI, etc. The influence of nonadia- 
batic matrix element, mixing excited and ground states, 
on the energies of resonances essentially depends on the 
phonon frequency. While in the adiabatic approximation 
ST transition is sudden and nonanalytic in A |23j |. nona- 
diabatic matrix elements turn it to smooth crossover |'26| . 
Thus, as illustrated in Fig. [SJi, the smaller the frequency 
the sharper the kink in the dependence of excited state 
energy on the interaction constant. 

Cautions should be made about approximate form of 
EPI J3J. Strictly speaking, actual momentum depen- 
dence of the interaction constant a (2^, can slightly 
change the obtained differences between nodal and antin- 
odal points though the general trends have to be left in- 
tact because ST is caused solely by the short range part of 
EPI [23| • Also the shell model calculation j2^| has shown 
that the dominant EPI is the short range one associated 
with the oxygen displacements. Therefore, even though 
the participation of Cu ions [2^] in the lattice vibrations 
relevant to EPI violates simple relation ~ 1/y/M, this 
relation is approximately valid. 

Finally, we discuss the relevance of the present results 
to the isotope experiment on the doped high T c cuprates 
• The most puzzling feature is that the IE is negligible 
near the phonon frequency, and is mainly observed in the 
high energy region. This can not be explained by the con- 
ventional Eliashberg-Migdal theory, which predicts the 
suppression of the multiphonon emission at high energy 
and the IE only near the one-phonon frequency. The in- 
termediate and strong coupling polaron theory gives a 
completely different picture, where the high energy part 
consists of the multiphonon sidebands and is subject to 



the IE as shown in this paper. Although it is not clear 
at the moment when the crossover from the doped hole 
picture to the metallic large Fermi surface occurs as the 
doping proceeds, the IE experiment Q suggests the for- 
mer might persists even at the optimal doping. However 
more study is needed both theoretically and experimen- 
tally for the doped case. In the undoped case, on the 
other hand, the present results can be directly compared 
with the experiments and more solid studies can be done. 
It is found that isotope effect on the ARPES lincshape of 
a single hole is anomalously enhanced in the intermediate 
coupling regime while can be described by simple inde- 
pendent oscillators model in the strong coupling regime. 
The shift of FCP top and change of the FCP height are 
relevant quantities to pursue experimentally in the inter- 
mediate coupling regime since IE on these characteristics 
is enhanced near the self trapping point. In contrast, shift 
of the leading edge is the relevant quantity in the strong 
coupling regime since this value increases with coupling 
as y/\. These conclusions, depending on the fact whether 
self trapping phenomenon is encountered in specific case, 
can be applied fully or par tially to another compounds 
with strong EPI [Tj, [H 13 • 

Note added in proof. Related studies of the isotope 
effect on ARPES were recently reported by S. Fratini 
and S. Ciuchi in the framework of dynamical mean 
field approach to the Mott-Hubbard insulator and by G. 
Seibold and M. Grilli |3lJ where coupling of quasiparticle 
to a charge collective mode is treated within the simple 
perturbative scheme. 
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